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Abstract 

Here we present a model to study the micro-plastic regime of a stress-strain curve. In this model 
an explicit dislocation population represents the mobile dislocation content and an internal shear- 
stress field represents a mean-field description of the immobile dislocation content. The mobile 
dislocations are constrained to a simple dipolar mat geometry and modelled via a dislocation 
dynamics algorithm, whilst the shear-stress field is chosen to be a sinusoidal function of distance 
along the mat direction. The sinusoidal function, defined by a periodic length and a shear-stress 
amplitude, is interpreted to represent a pre-existing micro-structure. These model parameters, 
along with the mobile dislocation density, are found to admit a diversity of micro-plastic behaviour 
involving intermittent plasticity in the form of a scale- free avalanche phenomenon, with an exponent 
and scaling-collapse for the strain burst magnitude distribution that is in agreement with mean-field 
theory and similar to that seen in experiment and more complex dislocation dynamics simulations. 
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I. INTRODUCTION 



In 1964, using state-of-the-art torsion experiments, Tinder and co-workers were able to 
achieve a strain resolution of 10 -8 for sub mm-sized samples to study the micro-plastic 
regime of highly pure poly-crystal Cu samples^, followed by tests on Zn single crystals some 
years later^. Their plastic-strain versus stress curves contained plateaus of stress which 
were attributed to the occurrence of discrete dislocation glide activity over a length scale 
comparable to the dislocation spacing of an assumed three-dimensional network. Indeed the 
authors write in refP, "The results suggest that an important fraction of the total strain, in 
the initial stages of deformation, involved motion of a few favourably situated dislocation 
segments through distances large enough to form new interactions with other elements of 
the three dimensional network. If this were so, then most elements of the network must 
have been relatively immobile, making little or no contribution to the strain." That most 
dislocations remain immobile remains a contemporary viewpointP. 

Another more recently pursued route is to probe discrete dislocation activity via the 
stress-strain curve of micron sized focused ion beam (FIB) milled single crystals. Here, 
advantage is taken of a nano- indent at ion platform equipped with a flat punch tip to com- 
press micro-/nano-crystals^^. In addition to the sub-nanometer displacement resolution of 
the system, which notably has a lower strain resolution than the above mentioned torsion 
experiments, the sample size is decreased to the micron range and below. As a result, the 
strain associated with the discrete dislocation activity is increased to an easily detectable 
magnitude. Whilst this more contemporary work has been primarily motivated by the 
"smaller-is-stronger" size-effect paradig m 1 4 1 5 1 8 ^ , an extensive analysis of the statistics of the 
discrete dislocation activity has revealed power-law behavior in the distribution of strain- 
burst magnitudes giving an exponent of ~ 1.6 — 2. 2 6 * 7 * 11 !. Similar exponents can also be 
found in bulk samples via detailed analysis of load displacement signals, where structural 
evolution is also seen to occult™ The very recent work of Dahmen and co-worker d 14 * 15 ! 
suggest that the variation in literature exponent values could be due to the different stress 
intervals used to bin the strain burst magnitude data, where a total integration over the 
stress interval to material value should, within mean field theory, give an exponent equal to 
precisely two. Such exponents are indicative of crackling^ or Barkhausen noise, and more 
generally of avalanche phenomena, indicating that dislocation mediated plastic deformation 
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belongs to a universality class that encompasses many natural phenomenon over a variety 
of different length and timescales. Indeed, similar power-law exponents can also be found 
for metallic glasses^ in which the underlying plastic deformation is fundamentally different 
to crystalline metals. 

Another class of experiments revealing the intermittent nature of dislocation dynamics 
is the acoustic emission monitoring of ic d 18 * 19 l Such experiments measure the acoustic en- 
ergy released by intermittent dislocation activity during constant stress deformation (in the 
tertiary creep regime). Indeed, via single sensor acoustic emission signals, such dislocation 
activity can be well characterised in time revealing power-law behaviour with exponents 
of 1.6 to 1.8, which is very similar to that seen via the micro-compression stress-strain 
curve ^6i7UU Furthermore, via multiple sensor monitoring, time and space clustering of dis- 
location avalanches could be observed^. It was found that avalanche epi-centres were corre- 
lated in space according to a non-integer power-law exponent indicating scale-free clustering 
and at short enough times such clusters were correlated in time indicating collective activity. 

That such scale-invariant dislocation activity occurs, is a signature of an underlying com- 
plex dislocation based micro-structure. An entity whose properties and evolution under 
an applied stress play a central role in the more general subjects of material strength and 
strain hardening 2 ^. Due to the complex dynamics and evolution of the dislocation structure, 
computer simulation based approaches have helped greatly over the past decades to clarify 
the underlying dislocation based mechanisms responsible for such structural evolution. One 
such method is the so-called dislocation dynamics (DD) approach. Early work involved two 
dimensional arrays of straight edge dislocations interacting via elasticity using single- and 
multi-slip geometries^ 221 . These works demonstrate that under an external stress, disloca- 
tion patterning emerges analogous to what is seen in static transmission electron microscopy 
experiments. Other and more recent works have developed these numerical techniques in 
terms of efficiency^!, strain boundary-conditions and obstacle/composite geometries^ to 
study more complex patterning such as the emergence of granular cell structures and the 
study of grain boundary network evolutional, strain hardening, and material strength in 
both bulk and confined volume systems^^. Analogous methodologies and investigations 
have also been carried out in three dimensions'^. 

Intermittent dislocation activity, in the form of dislocation avalanches, has also been 
studied using the DD simulation techniqu e! 18 * 28 * 31 ! Such simulations have produced power-law 
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exponents of the distribution of strain burst magnitudes similar to that seen in experiments, 
indicating some degree of scale-free behaviour, where dislocations are arranged in meta- 
stable cell and/or wall structures, with only a minor fraction of the dislocation population 
moving intermittently, thereby creating discrete strain jumps. The so-called self organised 
criticality (SOC) view of the dislocation network state offers one theoretical platform for the 
understanding of the observed universality^, in which the material system organises itself 
into a configuration that is critical, resulting in scale-free behaviour upon transiting to a 
new realisation of the critical state. Originally developed to describe sand-pile dynamics^, 
the approach is somewhat at odds with the historical viewpoint that dislocation structure 
evolution is primarily driven by equilibrium driving forces such as that embodied in Low 
Energy Structures (LES) theory^ and in a wide range of strain hardening theories. The 
finding that the occurrence of avalanche phenomenon is insensitive to the nature of the 
forming immobile dislocation network, the slip geometry, the deformation mode, and the 
details of the dislocation dynamics and spatial dimension is a central hallmark of SOC, 
which is robust against the details of the underlying physical model. 

This motivates investigating plastic flow with simpler models that explicitly do not take 
into account the fine details of individually interacting sessile and mobile dislocations, nat- 
urally shifting the focus of plastic flow from complex dislocation structural evolution to the 
interaction of a minute mobile dislocation population with a simplified description of the 
sessile dislocation population. Such an approach is analogous to the study of dislocations 
in the presence of pinning potential J 37 * 38 * and more generally to coarse grained models of 
plasticity that study the depinning transition (see ref.^ and references therein, and also 
refsJ 14 * 15 1 ). Dislocation field theories are also well able to study intermittency, however unlike 
dislocation dynamics based methods which are able to simulate only small plastic strains, 
these models are able to incorporate dislocation transport over non-neglible distances and 
therefore able to simulate significant structural evolution^. 

In this work, a simple model is therefore proposed in which a dipolar dislocation matP^ 
is embedded into an internal static sinusoidal stress field defined by a wave-length and 
a shear-stress amplitude. The explicit dislocation population is modelled via a standard 
DD algorithm. It is found that, similar to very complex and detailed 3D-DD simulations, 
the resulting stress-strain curves evolve in a discrete manner that reflects an underlying 
intermittent plasticity originating from irreversible changes in dislocation configuration, and 
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that the distribution of the corresponding strain-burst magnitudes reveal both extremal 
value statistics and scale-free avalanche behaviour. Thus, although the complex details of 
microscopic dislocation mechanisms and structure are omitted, the simple model is still able 
to capture the fundamental properties of intermittent micro-plastic flow. 

In the following section, the model is formally introduced and the DD technique described, 



and in sec. Ill various loading modes to produce a stress-strain curve are presented. Sec. IV 



presents the DD simulations for a wide range of model parameters to investigate their effect 
on deformation behaviour. The statistical analysis of the stress at which intermittent activity 
occurs, as well as the distribution of related strain-bursts, are presented in sec. [V) Finally, 
sec. 



VI discusses the context for the model, where it is argued that its applicable range is 



restricted to the micro-plastic region of the stress-strain curve — a regime where significant 
structural evolution and work hardening are, to a large extent, absent^ ^. 



II. DESCRIPTION OF MODEL 



The proposed model consists of two parallel slip planes, in the x — z plane, separated by a 
distance h along the y direction. Each slip plane is populated by infinitely long straight edge 
dislocations where one slip plane contains N + dislocations with Burgers vector b = (b, 0, 0) 
and the other N_ dislocations with Burgers vector b = (—6, 0, 0). Equal numbers of each type 
of dislocation are considered to ensure no net Burgers vector content: N + = N- = N/2. Such 
a structure is traditionally known as a simple dipolar matPl To this, a time independent 
internal sinusoidal shear-stress field is added of the form 

TlnternalO) = Tq COS (-^"J , (!) 

which is parametrised by a shear-stress amplitude To and a periodic length scale A. Fig. [T]dis- 
plays a schematic of the model system in which the line direction of each straight dislocation 
is perpendicular to the plane of the figure. The dislocation density is defined by the total 
number of dislocations divided by the area of the system. The area of the system is defined 
by its length along x, defined as d, and its spatial extent along y, defined as 2h. Detailed dis- 
cussion on the origin of the internal sinusoidal stress field is deferred to sec. |VI[ however, its 
existence should be viewed as a simplified representation of a pre-existing (and unchanging) 
immobile dislocation population. The explicit dislocations, schematically shown in fig. [TJ are 
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FIG. 1: Schematic of the two-dimensional dislocation dipolar mat system consisting of edge dislo- 
cations populating the two slip planes and the external sinusoidal stress field visualised in blue. 

therefore to be viewed as the mobile dislocation population of the system, and are quantified 
by their density p m . 

In addition to the force arising from the internal shear-stress field, each dislocation will 
experience a force of elastic origin from all other dislocations within the system. Since it is 
assumed that the dislocations are unable to leave their slip plane, the only relevant force per 
unit dislocation length component will be along the x direction and, via isotropic elasticity, 
may be calculated from the well known expression between two edge dislocations (labelled 
i and j) on parallel slip planes^ 

, (■'i'.J'r.j Ax(Ax 2 - Ay 2 ) 

Jx ' lj 2tt(1 - v) (Ax 2 + Ay 2 ) 2 ' 1 j 

Here G is the isotropic shear modulus and v Possion's ratio of the isotropic elastic medium, 
b x ,i (b x ,j) is the Burgers vector in the x-direction of the ith (jth) dislocation, and (Ax, Ay) is 
the two dimensional vector defining the dislocations' spatial separation. Presently a model 
isotropic Cu system is implemented, in which the shear modulus is taken as G = 42 GPa, 
the Possion's ratio as fi = 0.43, and the Burgers vector magnitude as b = 2.55 A. 

For the present work, periodic boundary conditions along d, the dipolar mat direction, and 
open boundary conditions along h, are assumed. Due to the long range nature of eqn. [2j the 
correct treatment of periodicity invovles the summation of all dislocation image contributions 
to the force per unit dislocation length on a given dislocation. For the considered one- 
dimensional periodicity, an exact solution to such a summation is tractable, and is given 
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Gb x ^b x j . f 2nAx 



by 

[d (cos (^) - cosh (g^)) + 2^Ay sin (^)] 
rf 2 (cos cosh (^)) 2 

A simple derivation of this equation is detailed in appendix [A] 

The temporal evolution of a particular dislocation configuration is characterised by the 
choice of an empirical mobility law. Due to the actual discreteness of the lattice at the 
atomic scale, a dislocation segment must overcome an energy barrier associated with the 
local shearing of atoms in order to move an atomic distance. This so-called Peierls energy 
barrier and the associated Peierls stress 4 ^, the stress at which the dislocation can begin to 
move (defined at a given temperature), results in the dislocation moving quasi-statically 
from atomic lattice site to atomic lattice site. At the meso-scopic scale this results in over- 
damped motion where the dislocation's velocity is proportional to the force acting on the 
dislocation — the present mobility law. The material specific proportionality constant is 
referred to as a damping parameter and is dependent on dislocation type, geometry, and on 
temperature. 

The equation of motion along the x direction for the ith dislocation is then given by 

5xi F, 



X,l 



u b ■ < 4 > 

where B is the damping coefficient, which for Cu is 5 x 10~ 5 Pa ^ 23 | 44 [ In eqn. |4j F Xj i, is the 
total force per unit dislocation length acting on the dislocation, 

F X) i [^Internal ^External] ~\~ ^ fx,ij- (5) 

Here TExtemai is an externally applied homogeneous shear-stress field and ri ntem ai is the static 
sinusoidal shear-stress field defined in eqn. [Tj 

The numerical solution of eqn. [4] constitutes the Dislocation Dynamics (DD) algorithm 
presently used in which an appropriate finite time-step, 5t, is used to integrate the equations 
of motion. The corresponding shear-strain response 5e within this St is calculated via 

1 N 

fc =^5>^- ( 6 ) 
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It is again emphasised that over the time-scale of St all atomic scale aspects are averaged over 
and inertial effects are ignored. Since the edge dislocations are infinitely long and straight 
such dynamics falls into the class of two dimensional DD modelling. 

The dislocation density is given by p m = Mr. For the present work, the periodic length 
d is chosen to define the distance h between the two populated slip systems via 2d/N. This 
sets the mean distance between dislocations along the x— and y— direction to be the same, 
and defines the dislocation density as (N/2d) 2 . Thus, h = and choosing a value of 

the dislocation density will fix the scale of interaction between the two parallel slip planes of 
the dipolar mat geometry — see appendix A. The motivation for such a restriction is to give 
the dislocation density a greater bulk-like relevance, where in the bulk limit it represents an 
isotropic density of dislocations. Other definitions of h are, of course, possible. 

III. LOADING AND THE CALCULATION OF STRESS-STRAIN CURVES 

To simulate a loading experiment and thus a stress-strain curve, a sample must be pro- 
duced and a loading mode chosen. Presently, sample preparation involves the chosen number 
of dislocations being initially placed at random positions within the dipolar mat geometry, 
and the structure relaxed via eqn. [4] to minimise the force on each dislocation to within 
a chosen tolerance. This is performed at r Externa i = 0. In the present work, no attempt is 
therefore made to determine a low (or lowest) energy dislocation structure, an approach that 
is compatible with the fact that the explicit dislocations within the model are only those of 
the mobile variety. This aspect is found to be crucial to the properties of the model and will 
be discussed in more detail in sec. [VTl 

There exist a number of ways in which a deformation simulation can be done. The first 
such loading mode is referred to as "stress-relaxed" and involves incrementing the external 
shear-stress by a value St (rExtcmai —> ^External + St), and relaxing the structure until the the 
sum of the dislocation force magnitudes, or equivalently, the dislocation velocity magnitudes, 
varies by less than the fraction 1CT 8 . 

During the relaxation, the associated plastic shear-strain increment may be calculated as 
the sum of eqn. [6] over all St time steps of the relaxation. Once convergence is obtained this 
shear-strain increment is added to the total plastic shear-strain and the cycle is repeated 
until the desired stress-strain curve is obtained. Fig. [2^i displays a typical shear-stress versus 
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plastic shear-strain curve during the initial stages of loading. Discrete dislocation activity, 
via strain bursts, is clearly evident and is separated by continuous an-elastic regions. By 
adding the elastic shear-strain, r/G, to the plastic shear-strain, the total shear-strain is 
obtained. Fig. [2]d displays the corresponding shear-stress versus total shear-strain curve for 
a small range of stresses and an approximately linear stress-strain curve with intermittent 
stress plateaus is seen. When displaying a similar curve for the full range of stresses, as seen 
in fig. [2^,, only a straight line is resolvable indicating that such strain bursts are well within 
the micro-plastic regime of deformation. At larger stresses a plastic flow regime is entered 
which will be investigated in more detail in subsequent sections. 




plastic shear strain total shear strain 



FIG. 2: a) A typical stress versus plastic strain curve produced using the "stress relaxed" loading 
mode. Stress plateaus indicating strain burst activity is clearly evident, b) Displays, for the same 
simulation, the corresponding stress versus total (elastic plus plastic) strain for a small range of 
shear-stresses. 

Experimentally, two distinct deformation modes can be used: displacement controlled 
and load controlled. Here we consider an inherently force controlled testing device. In such 
a case, displacement controlled testing is done by adjusting the applied load via a feed- 
back loop such that the displacement rate is held at a fixed value throughout the loading, 
whereas for load controlled experiments, the applied load simply increases at a chosen rate. 
For the present model, these deformation modes correspond to a constant shear-strain rate 
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and constant shear-stress rate loading condition. To obtain a stress-strain curve with a 
constant shear-stress rate, a numerical value for the applied stress rate, r Externa i, is chosen. 
This then defines a stress increment St = f ExtCTnSL iSt, where St is the time-step used to evolve 
the dislocation network according to eqn. |U Thus at every simulation iteration, the stress 
is increased by St and the configuration evolves in time by an amount St. To implement a 
constant shear-strain rate loading mode, e, a numerical value is chosen and the appropriate 
St stress increment, to achieve such a strain rate, is performed every simulation step. The 
actual value of St is determined by assuming that the total strain rate decomposes additively 
into an elastic and plastic component: 

. . . _ ^"External . /-a 

£ ^elastic ~r Aplastic ^ t Aplastic; \' ) 

giving 

^External = G (i — Aplastic) (8) 

or 

St = G (e - plastic) St, (9) 

where e p i a stic^ is the strain increment per simulation iteration, calculated via eqn. [6j Thus 
the stress increment St is determined by the correction needed to achieve the required 
constant strain rate for the next simulation iteration. In the above, due to the simplified 
geometry of the model (fig. [I]), a pure shear modulus, G, rather than a Young's modulus is 
used. 

To investigate how such loading modes affect the discrete dislocation activity seen in 
fig. [2j the appropriate stress-rate and strain-rate must be chosen. This is done by first 
choosing e, from which ^External is obtained via e/G to ensure that in the elastic/an-elastic 
regime both deformation modes have the same total strain rate. Fig. [3] displays a single 
strain burst in all three considered loading modes. For the constant strain rate and stress 
rate modes two strain rates are considered: 0.1 s _1 and 1.0 s _1 . For the "stress-relaxed" 
loading mode a sharp plateau is evident in fig. |3^ with an identically zero gradient. In this 
region, the constant stress rate mode also exhibits a plateau but with a (non-zero) positive 
gradient since, during the evolution of the strain burst, the stress is rising at the chosen rate. 
Also, the stress at which the strain burst initiates is somewhat higher (and increasing with 
increasing strain rate) than that in the stress relaxed mode indicating a strain-rate effect. In 



10 



a) 



C3 



24.64 — 



24.62 — 
S3 24.6 
24.58 h- 



X 

CO 




1 1 I r 

stress relaxed 
stress rate (1.0 s L ) 
stress rate (0.1s ) 
strain rate (1.0 s ) 
strain rate (0.1s ) 



x x^ x x^ x x^ xx^ 

plastic strain (%) 



24.56 



5 V V 6 V V 6 V V 6 
total strain (%) 



FIG. 3: Close up of a dislocation burst for the three considered loading modes and at different 
strain rates, a) Stress versus plastic strain and b) stress versus total (elastic plus plastic) strain. 

this regard, the "stress-relaxed" deformation mode can be considered as the zero stress-rate 
limit of the constant stress rate loading mode in which the dislocation configuration always 
has time to relax before the next stress increment. For the constant strain-rate loading 
mode, the onset of the strain burst occurs at similar stresses to that of the constant stress 
rate, however as the strain-burst evolves, the stress decreases to maintain the chosen strain- 
rate. Fig. [3)3 displays the same strain burst with the stress now as a function of the total 
strain. The greatest effect is seen in the constant strain-rate mode where due to the drop 
in stress during the strain-burst, there is a rapid drop in elastic strain. Fig. displays the 
plastic strain rate as a function of plastic strain for the burst shown in fig. [3] (for a constant 
strain rate of 1.0 s _1 ). The elastic strain rate response as dictated by eqn. [9] is also shown 
and correspondingly reduces to compensate for the rise in the plastic strain rate. Data for 
a constant strain rate of 0.1 s _1 differs little from the 1.0 s _1 data, a result (along with the 
high strain-rates of fig. [4^,) due to the driven zero temperature nature of the simulation. 

It is noted that for the constant strain rate loading mode, the loading system responds 
instantaneously (to within St) to any discrete plastic event (hence leading to the backwards 
curvature seen in fig. ^p)- To model an instrumentally realistic device, a more complex 
differential equation than that of eqn. M would need to be developed, which takes into 
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account the finite delay time and resolution of the loading apparatus. 




1.12e-05 1.14e-05 1.16e-05 1.18e-05 5000 10000 15000 



plastic shear strain simulation iteration 

FIG. 4: a) Elastic and plastic strain rates as a function of strain for the dislocation burst displayed 
in fig. [3j b) Average dislocation velocity during the same dislocation burst in "stress relaxed" 
deformation mode. 

Fig. |4]d now shows the average dislocation velocity during the "stress relaxed" defor- 
mation mode for four consecutive applied shear stresses in which the first constant stress 
relaxation undergoes the irreversible plastic strain rate seen in fig. |3j For this stress value, an 
initial relaxation of the dislocation velocities is seen, followed by a rapid acceleration of the 
dislocations corresponding to the emerging instability and the eventual irreversible plastic 
event of fig. [3] This continues until a maximum velocity is reached at which the dislocations 
then begin to decelerate until the convergence criterion is met. Three further constant stress 
increments result in an immediate relaxation of the dislocation velocities corresponding to 
a reversible (anelastic) relaxation of the dislocation configuration. 

IV. RESULTS 

In this section the "stress-relaxed" loading mode is used to investigate the influence of 
the model parameters on the stress-strain curve. 



12 



A. Stress-strain behaviour as a function of micro-structural parameters: To, A, 
and the mobile dislocation density, p m . 



The influence of the sinusoidal shear-stress field amplitude is first investigated. Fig. [5^ 
displays the resulting shear-stress versus plastic shear-strain curves for three different choices 
of To; namely lOOMPa, lOMPa and IMPa for a system with d = 200 fim and A = 2 /im, 
containing N = N + + N_ = 20 + 20 dislocations. This gives a mobile dislocation density 
of p m = 1 x 10 10 m -2 . The figure demonstrates that the choice of r strongly controls the 
stress at which macroscopic plastic flow occurs. In fig. [5^i, the vertical axis is plotted as a 
logarithmic scale to reveal the fine structure of the micro-plastic region. In all cases discrete 
strain bursts are evident. Fig. [5]d displays the corresponding shear-stress versus total shear- 
strain curves. The small plastic shear-strain values evident in figure [5^i and the sharp yield 
transition in [5]d emphasise that the shear-stress versus plastic shear-strain data prior to flow 
is clearly in the micro-plastic regime of the stress- strain curve. 
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FIG. 5: a) Stress versus plastic strain for three different To values, for a mobile dislocation density 
1 x 10 10 m -2 , and b) the corresponding stress versus total strain. 



In figs. [6^i-b, p m is now varied from 1 x 10 10 m~ 2 to 2.56 x 10 12 m -2 by changing the 
number of dislocations for a fixed d = 100 /xm and A = 2 //m. Here To is 50 MPa. Figs. |6^i-b 
demonstrate that with increasing mobile dislocation density the flow stress decreases from 
tq. Thus at low enough dislocation densities (as in fig. [5]), t defines the flow stress of 
the system, a natural result since a well isolated dislocation will be primarily affected by 
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the sinusoidal shear-stress field. Figs. also demonstrate that the micro-plastic regime 
broadens with increasing dislocation density and that the first discrete dislocation event 
occurs at a decreasing applied homogeneous shear-stress (see arrowed regions in fig. |6|d). 
Figs. [6}>d show similar data for tq = 10 MPa, for the same values of dislocation density. 
Although a similar trend is seen, the deviation away from the maximum flow stress of tq = 10 
is a proportionally stronger function of the increasing dislocation density. 
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FIG. 6: Stress versus plastic strain for different mobile dislocation densities for (upper panels) 
To = 50 MPa and (lower panels) tq = 10 MPa. The right panels display the same data with strain 
plotted using a logarithmic scale to emphasise the initial micro-plastic strain region where the 
arrows (in b) indicate the first occurrence of a strain response. 
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The remaining micro-structural variable that can be varied is the characteristic length 
scale of the sinusoidal stress field, A. Fig. [7] displays stress versus plastic strain curves as A 
increases from A = 0.5/im to A = 20/im, with r = 10 MPa and 50 MPa, d = 800/im and 
p m = 1 x 10 10 m~ 2 . Here a larger d was used to obtain a better statistical sample of the 
dislocation environment at the length scale of A. The chosen dislocation density is similar 
to that of fig. [5] in which tq largely determined the flow stress. Inspection of fig. [7] reveals 
that with increasing A the micro-plastic regime is broadened due to increasing strain burst 
magnitudes with increasing A. For small strain (figs. [7^l and up), it also appears that the 
flow stress regime occurs at a reduced stress (for increasing A), however at larger strains 
(figs. [TJd and [7]l) flow stresses approximately equal to r are eventually reached. It is also 
noted that with increasing A the strain burst magnitude increases. With decreasing tq, the 
deviation of the curves away from To is a stronger function of increasing A resulting in a 
somewhat broader micro-plastic region. 

The central properties of the present model are reflected in figures [5j [6] and [7} If the 
dislocation density is low enough or r sufficiently high, then r will primarily determine 
the shear-stress at which the extended plastic flow regime begins. In this limit, the internal 
shear-stress that each dislocation experiences is due to the sinusoidal stress field and when 
^applied approaches r , the plastic flow regime is immediately encountered with a negligible 
micro-plastic region. By increasing the mobile dislocation density or decreasing To the micro- 
plastic regime broadens and also the stress at which extended plastic flow occurs decreases 
due to the increasing influence of the elastic interaction between the dislocations. This latter 
feature is also reflected in the decreasing stress at which the first strain burst is seen. By 
increasing A, the strain burst magnitude increases which broadens the micro-plastic regime. 

B. Stain burst behaviour as a function of periodic length scale 

The stress versus plastic strain behaviour for three different values of d, using To = 10 
MPa and p m equal to 1 x 10 10 m~ 2 is presently investigated. For all three simulations A = 2 
fim. Fig. [8] displays the resulting stress-strain curves demonstrating that there is no strong 
overall d dependence. Indeed, the stress at which the first strain burst occurs differs little 
for the three samples: 3.4 MPa for d = 400/iin, 3.2 MPa for d = 800/im and 3.0 MPa for 
d = 1600/im. It is seen, however, that with increasing d the scale of the discrete plastic 
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FIG. 7: Stress versus plastic strain for a range of A choices, where for each simulation tq = 50 
MPa (upper panels) and tq = 10 MPa (lower panels), d = 800/um and a mobile dislocation density 
1 x 10 10 m~ 2 . The left panels show the initial strain regime and the right panels show a larger 
strain range. 

strain bursts becomes finer, leading to small-scale differences in the curves. Inspection of 
the minimum magnitude of the plastic strain bursts observed in this figure, and in figures [5] 
and|6j reveals it to be approximately |6| X/2dh, indicating that the strain bursts correspond 
to dislocation motion over a distance equal to at least A. It is worth noting that the existence 
of such minimum discrete strain bursts is an artifact of the periodicity used to approximate 
a system of infinite extent. Indeed, in the limit d — > oo the minimum strain bursts would 
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approach zero. 
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FIG. 8: a) Early stress versus plastic strain for different periodicity lengths d with To = 10 MPa 
and a mobile dislocation density 1 x 10 10 m -2 . b) A close up of the stress-stain curves indicating 
that with increasing periodicity the minimum discrete strain bursts reduce in strain magnitude. 

V. STATISTICAL PROPERTIES OF INTERMITTENT FLOW 
A. First burst behaviour 

To better understand the relationship between the mobile dislocation density and A, a 
series of simulations is performed to investigate the statistics of the shear-stress at which 
the first strain burst occurs. For computational efficiency these simulations were done using 
the constant stress-rate loading mode with d = 800 //m and r = 50 MPa. Fig. [9^i displays 
such shear-stress values as a function of A for three different mobile dislocation densities. 
For each A and dislocation density value, the results of many different simulations are shown 
to indicate the degree of statistical variation. Inspection reveals that by increasing the 
dislocation density, the first-burst stress value becomes increasingly dependent on the value 
of A, decreasing with decreasing A. For all three considered dislocation densities, the degree 
of scatter increases with decreasing A. 

The state of the dislocation configuration prior to loading provides an understanding 
of the results of fig. [9^,: N dislocations are introduced into the dipolar mat by randomly 
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FIG. 9: a) Plot of stress of first strain burst as a function of A for three dislocation densities for 
T = 50 MPa and b) r = 10 MPa at a dislocation density of 0.1563 x 10 10 m- 2 . 



placing them along the two possible slip planes and the system is then relaxed to a local 
equilibrium configuration. For large enough A this relaxation process will result in similar 
numbers of dislocations for each A unit, producing a dislocation configuration that is in- 
creasingly ordered. However, as A decreases, there will exist an increasingly varying amount 
of dislocations for each A unit, and the original disorder associated with the random initial 
positions of the dislocations is increasingly preserved. These two trends constitute the origin 
of the increased scatter and general decrease seen in fig. [9^ with respect to decreasing A, 
since the first-burst stress value is directly related to the lowest critical shear-stress required 
for an irreversible configurational change to occur. That is, the first-burst stress probes 
the extremal values of the dislocation environment and a greater variation in dislocation 
environment will naturally lead to a greater variation in the first-burst stress, which also 
increases scatter and results in an on-average reduced stress magnitude. Dislocation inter- 
action enhances this effect by generally increasing the degree of disorder in the low-lambda 
limit. 

The above result is, in fact, dependent on the system size (the periodicity length, d) being 
finite. As the system size grows at a fixed value of A and dislocation density, the chances 
for an extremal dislocation configuration to occur increases, resulting in a increased scatter 
and reduction of the first-burst shear-stress. It is therefore expected that as the system size 
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increases, the saturation in first burst stresses as a function of increasing A will be shifted 
to larger values of A. Such a size effect must also be reflected in the choice of r which, 
along with the dislocation interaction, determines the degree of configurational disorder in 
the initial state. For a given value of A, a smaller value of tq should increase the disorder due 
to the increasing influence of the dislocation interaction on the initial configuration. Fig. 
displays a similar figure of first burst stresses as a function A using a value of tq = 10 MPa 
and, indeed, demonstrates a shift to the right of the saturation region. 

Although these results are an artifact of the system size, they reveal that through proper 
choice of model parameters, the intermittent plasticity seen in the past sections is driven 
by the extreme value statistics associated with the initial dislocation configuration — a 
behaviour that is to be equally expected in real crystals. 



B. Distribution of strain burst magnitudes 

The statistical properties of the strain burst magnitudes are now investigated for a system 
with To = 10 MPa, A = 2 /im, and p m equal to 1 x 10 10 m~ 2 . Three systems with d = 400, 800 
and 1600 /im are first considered. The stress-strain curves for the d = 400/im and d = 800/im 
systems span up to the flow stress regime at approximately r = 10 MPa, whereas for the 
d = 1600/zm data is only available up to approximately 9.4 MPa. Fig.[l0]plots the magnitude 
of the strain burst as a function of the shear-stress at which it occurred, for all three samples. 
In the first instance, the plastic strain magnitudes for each stress step (as calculated via eqn. 
6) were analysed and it was found that irreversible plasticity via discrete strain bursts begins 
at approximately 3 MPa with magnitudes equal to |6| X/2dh (indicated by the corresponding 
coloured arrows). With increasing stress, larger strain bursts occur with magnitudes equal 
to multiples of the basic strain burst unit for each sample. Fig. [TUJd displays the entire range 
of strain bursts using a vertical log-scale, and it is seen that the spread in strain bursts 
increases rapidly as the shear-stress approaches the flow regime at approximately tq = 10 
MPa. 



Plotting the strain-burst magnitude distribution derived from the data of fig. 10, reveals 
a power law behaviour which is similar for each of the three considered systems. However 
from these graphs (not shown), a reliable exponent could not be obtained and to investigate 



this aspect further improved statistics is required. Fig. Hi shows strain-burst magnitude 
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FIG. 10: a) Magnitude of strain bursts as a function of the applied shear-stress for three samples 
differing only in their size d. The arrows indicate the (minimum) fundamental strain bursts for 
each sample, indicating that larger strain bursts occur as multiples of these values, b) Vertical log 
plot for the entire applied stress range demonstrating that as the flow stress regime is approached 
individual strain bursts can become large. 

distributions from a much larger number of events, plotted using a log-log scale, for three 



different choices of parameters: 1) A = 2 //m and tq = 10 MPa (same as that of fig. 10), 2) 
A = 10 /im and tq = 10 MPa, and 3) A = 2 fim and tq = 50 MPa. For all simulations d = 800 
fim and p m = 1 x 10 10 m 2 . For 1), eleven deformation simulations were performed, and for 
2) and 3) sixteen independent simulations were performed for each system. All distribu- 
tions were obtained by logarithmic binning of identified irreversible plastic strain events. 
Inspection of this figure reveals good power law behaviour over a number of decades. This 
is particularly the case for the parameter set 1), in which more statistics could be obtained 
from each individual curve. Indeed, for this parameter-set, fig. [9] shows that the first-burst 
stresses are low and exhibit strong scatter indicating significant intermittent plasticity. For 
this parameter set, a power law fit with an exponent of —1.96 ± 0.03 was found. Visual 
inspection of the remaining two parameter sets reveals an approximately similar exponent. 
For parameter set 2), finite size effects for the largest strain burst magnitudes are apparent. 
This is to be expected since at A = 10 /zm, there exists just 80 A lengths within the d = 800 
/j,m system, whereas for parameter sets 1) and 3), where A = 2, there exist 400 A lengths. 



The above mentioned data of fig. 11 1 was obtained by logarithmic binning of all discrete 
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irreversible plastic events in the micro-plastic regime. In the recent mean-field theory (MFT) 
work of Dahmen and co-worker d 14 * 15 ! this is referred to as the stress integrated distribution, 
Ant (SO in which S is the strain burst magnitude, whose integrand is given by the stress 
dependent distribution function, 

D(r,S)^S K f{S(r c -r) 1 ^), (10) 

where r c is the critical stress at which extended plastic flow occurs, and the exponents k 
and a, and the scaling function / are universal. From MFT theory k = 3/2 and a = 1/2. 
Performing the stress integral gives Ant (S) ~ S~^ +a) = S~ 2 where the last step uses the 



MFT exponents. The data of fig. Hi is therefore in excellent agreement with that of MFT. 



Further insight can be given via the stress-binned complementary cumulative distribution, 

POO 

C{t,S) = / dSD(r,S), (11) 



whose stress integral gives a power law ~ S^ K+a ^ = S 1 . Fig. 11 d, plots this quantity 



for the case of A = 2 /im and r = 10 MPa and gives an exponent equal to —0.99 ± 0.02. 



Fig. 11 3 also plots C(t,S) for a selection of stress intervals revealing strong cut-off effects 
due to the tunable criticality of the MFT model. Following refP^, these curves are found 
to collapse when appropriately scaled via the function / = (r c — r)/r c onto the universal 



mean field result, as shown in fig. [Tip , giving quantitative evidence that scale free avalanche 
behavoir does indeed occur for the current model. Here r c is taken as r . Similar behaviour 
is seen for the case of A = 10 /im and r = 10 but with / = (r c — t)/t c — c in which the 
c = 0.19 is an adjustment parameter to correct for finite size effects 15 . For the case of A = 2 
/jm and To = 50 the statistics was not of a sufficient quality to perform the procedure. 

When compared to more realistic two and three dimensional simulations that include 
dislocation reactions and much larger dislocation numbers, it might seem remarkable that 
such a simple model is able to exhibit scale-free behaviour. That this is the case, is a central 
feature of the SOC phenomenon in which the universal features of the critical dislocation 
configuration are robust against the details of the underlying physical model. Moreover, due 
to this simplicity, extremely good power law behaviour is obtained over many decades for a 
minimal computational effort. 
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FIG. 11: a) Statistics of strain bursts derived from stress-strain simulations for three choices of 
parameters. Data is shown as a log-log plot of the histogram of strain bursts in units of the 
fundamental strain burst magnitude, b) Corresponding stress dependent and stress integrated 
complementary cumulative distribution for the case of A = 2 fim and To = 10 MPa and c) via an 
appropriate scaling procedure where / = (tq — t)/tq all stress dependent complementary cumulative 
distributions collapse on the universal curve given by mean field theory — see refP^l. 



VI. DISCUSSION 



The simulations of sec. II VI demonstrate that two dominant factors control the characteris- 
tic stress scale of the stress-strain curve. These are the choice of To and the mobile dislocation 
density. Figs. [5] and [6] demonstrate that t sets an upper stress limit for extended plastic flow 
to occur. This upper limit is approached when the mobile dislocation density is sufficiently 
low that the sinusoidal stress field is the dominant contribution to the stress field each dis- 
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location feels. By increasing the mobile dislocation content, the characteristic shear-stress 
scale reduces due to the increasing role of the internal stress fields arising from the elastic 
interaction between dislocations. At the same time a broadening of the micro-plastic regime 
is observed, which very much is a general property of micro-yielding, where a greater num- 
ber of mobile edge dislocations correlates with larger measurable plastic strain^. Figs. [7] to 
9] demonstrate the effect of A (and d) on the stress-strain curve is somewhat subtler. The 
choice of A influences the initial configuration of the mobile dislocation population, thereby 
affecting the statistics of the stress at which the first strain-burst occurs and also the way in 



which the extended plastic flow regime is reached. Fig. 11 shows, however, that the statistics 
of the strain burst magnitude is insensitive to all three discussed parameters, reflecting the 
universality of SOC with respect to micro structural details. 

How shall one quantitatively interpret the imposed sinusoidal stress field? In a real mate- 
rial that is nominally free of dislocation structures, some type of length scale will naturally 
emerge as a function of macroscopic plastic strain due to the evolution of a growing and 
interacting dislocation structure that eventually leads to the phenomenon of patterning^!. 
Although patterning is a term predominantly referring to the effects of latter stage II and 
III hardening regimes, the emergence of micro-structure length scales is expected to occur 
at all stages of plasticity ranging from slip, dipole and eventually cellular pattern^ 36 * 47 * 48 !. In 
fact, a micro-structural length-scale can equally well be defined for an undeformed as-grown 
material, where the mean dislocation spacing can be used to describe the initially present 
internal stress fluctuations — a view point that is central to the early work of Tinder and 
co-workers^l. From this perspective the imposed sinusoidal stress field, can be viewed as 
the simplest realisation of the internal stress field arising from such a structure. 

In the present model this stress field is time independent, implying it is constructed 
by that part of the dislocation population that is immobile, with the explicit dislocations 
and their dynamics, arising from the (much smaller) mobile component of the dislocation 
population. Thus a typical loading simulation can be seen as the deformation of a model 
material that has a particular sample preparation or deformation history characterised by 
To and A, and a mobile dislocation density that is only a small part of the total dislocation 
density. 

Much past work exists concerning the emergence of internal stress and length scales 
as a function of deformation history. In early work on the theory of cell formation, two 



23 



relationships have emerged in which the total evolving dislocation density, ptotai, plays a 
central role. They are: 

''"material 7 , /i o\ 

q OC OVAotal (12) 

and 

Amateriai OC , (13) 

V Ptotal 

where G is a representative (not necessarily pure) shear modulus and b the Burgers vector 
magnitude. In the above, r mate riai is the evolving flow stress of the material and A mate riai is an 
evolving internal length scale that can be referred to as a cell size. The first expression has its 
earliest origins in a Taylor hardening picture in which the total dislocation density is seen as 
an immobile forest dislocation population. The second equation has its theoretical origins 
in the early work of HollP^ who derived it for a dipolar population of screw-dislocations, 
showing that a uniform arrangement was unstable to fluctuations with one such length scale 



dominating, characterised by eqn. 13 This length scale, which could be related to a fixed 
self-screening distance of the dislocation network, was postulated to reflect an emerging 
cell size. The approach was based on an energy minimisation principle, however due to 
dislocation reactions, the more modern viewpoint is that the dynamics of cell formation lies 
in a statistical process involving dislocation reactions and that the screening length, and 
therefore cell-size, is an evolving variable®'. 

Thus, the model parameter A has a direct counterpart in cell formation theory, A mate riai> 
which can represent quite generally, a mean-free path for a mobile dislocation, a dipolar 
screening length or a well evolved cell length scale. Moreover, since an unloading/loading 
cycle will generally return a system to the flow stress before unloading, and that the present 
simulations have shown that the flow stress is partly controlled by r , r should be in some 
way related to r materia i. From this perspective, r and A are parameters that are not entirely 



independent from each other. In fact, eqns. 12 and 13 express that the cell size decreases 
inversely as a function of flow stress, a well known experimental observation that is refered 
to as "similitude" 1311 . 

Although similitudity is generally confirmed by experiment, some experimental work does 
present a more complicated picture. Early tensile/TEM work on tapered Cu single crystals 
finds an initially broad distribution of cell sizes that narrows and shifts to small lengths 
with increasing flow stress^. This result suggests that a single structural length-scale might 
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not always be a good statistical description of the evolving micro-structure. Indeed, more 
modern viewpoints, in which dislocation structure evolution is a non-equilibrium process 
tend to suggest a distribution of emerging length-scales leading to a scale-free fractal-like 
structure. Although such micro-structures have been quantitatively established by TEM 
investigations of latter-stage hardened single crystals of Cu | 52 | 53 | ; their existence is not uni- 
versal, depending strongly on material type and deformation history. The current work does 
not address this aspect. More general forms of an inhomogeneous internal stress field that 
capture such scale-free micro-structural features can be envisioned; a direction which will 
be investigated in future work. 

Whilst the dipolar mat geometry in an external field offers a platform with which to study 
the depinning transition and more generally the transition to extended plastic flow, when 
comparing to experiment, careful consideration has to be given to its regime of applicability. 



To do this, a typical simulation of sees. IV and [V] is now broadly summarised. Upon choosing 
numerical values for all model parameters, the N dislocations are introduced to the system 
via a distribution of random positions. This unstable configuration is then relaxed to a local 
minimum energy in which the forces on each dislocation are below a small threshold value. 



The deformation simulation is then begun using one of the three loading modes of sec. Ill 



As the stress increases, intermittent plasticity increasingly occurs until a stress is reached 
at which extended and overlapping strain events occur, which in the previous sections has 
loosely been referred to as the plastic flow regime. 

It is important to emphasise that no attempt has been made to obtain the global energy 
minimum of the starting configuration. Such an initial state turns out to play a crucial role 
in the observed properties of the model, since many high energy configurations will exist, and 
it is these that dominate the early stages of plasticity. As a deformation simulation proceeds, 
such high energy configurations structurally transform eventually leading to a plastic flow 
regime and often to the homogenisation of the dislocation configuration. In other words, 
the extended plastic flow regime should be considered to be outside the applicability regime 
of the present model when a comparison to experiment is made, or equivalently, the present 
model is only suitable for the study of the micro plastic regime of the stress-strain curve. 

The rational behind the use of an initial high-energy dislocation configuration originates 
from the assumption that the explicit dislocation population of the model represents only the 
mobile dislocation network, which constitutes only a small part of the true population. Thus, 
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in the same way as To and A characterise the sample preparation or deformation history of 
the model material, so does the initial high energy (explicit) mobile dislocation content. This 
is quite compatible from the perspective of SOC in which the dislocation structure reaches 
a critical configuration that is far from equilibrium, and that structural rearrangements 
correspond to the system transforming from one SOC state to another. By construction, that 
part mediating the structural transformation will be the current mobile dislocation content. 
The central simplification of the present model, is that it separates the mobile and immobile 
populations, associating the former to an explicit mobile dislocation content that represents 
the non-equilibrium component of the network, and relegating the latter to an effective 
static internal stress field. That this internal stress field is unchanging and that the same 
explicit mobile dislocation population exists as a function of strain for the entire deformation 
simulation, is of course different from a real material, where the structure evolves with 
strain, and at any particular non-negligible strain interval, quite different dislocations might 
constitute the mobile dislocation population. This again emphasises that the present model 
should only be applied to the micro plastic regime, where significant structural evolution is 
minimal. 

Experimental evidence for a lack of structural evolution in the mico-plastic regime is best 
seen in low amplitude cyclic deformation experiments of FCC metals, in which the plastic 
strain per cycle can be as low as ~ 1CT 5 leading to significant changes in load stress and 
internal length scale only after the occurrence of several tens-to-hundreds of thousands of 
cycle^SES j s further noted, that documented experimental studies on the micro-plasticity 
at room temperature primarily report on movements of edge or non-screw type dislocations, 
whereas a clear increase in dislocation density or the formation of dislocation structures as 
a result of multiplication remains absent™ 153 . In the bulk case there are exceptions to 
this trend where in the case of a work-hardened Al-Mg alloy which exhibits dynamic strain 
ageing, emerging structural length-scales were already detected in the micro-plastic regime^ 
using high-resolution extensometry methodj 12 * 57 * 58 l 



The results of sec. VB demonstrate that the developed model exhibits power-law be- 
haviour in the distribution of strain burst magnitudes, and thus the scale-free avalanche 
phenomenon seen in experiments, either via the stress-strain curve of micro-compression 
tests^ or via in situ acoustic emission experiment J 18 * 19 *, and in simulation, via two or three 
dimensional dislocation dynamics simulations in which the entire network is represented by 
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an explicit dislocation population and individual dislocation reaction mechanisms are taken 
into account. With an exponent of approximately -1.96, the model gives a value that is 
somewhat higher than that seen in simple metals 6 and ice 18 , but more comparable to that 
seen in LiF crystals^. That such a simple model can admit scale-free behaviour, is connected 
to the dependence of the intermittent plasticity on the extremal configurations of the explicit 
dislocation population. This was directly seen in the statistics of the first-burst shear-stress 
and also the distribution of strain burst magnitudes, where with a large enough increase 
of A (say from 2 /im to 10 /im) the first burst statistics changes from being dominated 
by extreme value statistics to that being dominated by the statistics of the most probable 
(fig. ^p) corresponding to an increased presence of cut-off effects in the statistics of strain 
burst magnitudes (fig. [TT^i). This is a natural result of the observation that quantities that 
depend on extreme value statistics can exhibit power-law behaviour in their distributions, 
emphasising a connection with SOC that is related to only the mobile dislocation popula- 
tion being in a non-equilibrium state and not to the characteristics of the present simplified 
immobile dislocation network — a manifestation of a scenario referred to as "nearly critical" 
or "robust critical"^. 



VII. CONCLUDING REMARKS 



A simplified two dimensional dislocation modelling framework has been introduced in 
which the explicit interacting dislocation population, constrained to a simple dipolar mat 
geometry, represents only the mobile dislocation density component of the total dislocation 
density, and the much larger immobile dislocation population is described by a static internal 
sinusoidal shear-stress field defined by an internal shear-stress amplitude and wavelength. 
These model parameters, along with the initial non-equilibrium explicit mobile dislocation 
content characterise either the deformation or sample preparation history of the model ma- 
terial. Because of the static nature of the internal field and the lack of dislocation-dislocation 
reactions, upon loading, the present model is restricted to the micro-plastic region of the 
stress-strain curve, and therefore to a deformation regime for a given material that involves 
negligible structural evolution. Despite the simplicity of the model and the restriction to 
the micro-plastic regime, the deformation behaviour exhibits a rich variety of properties as 
a function of the model parameters. In particular, intermittent plasticity is observed whose 
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strain burst magnitude distribution exhibits scale-free avalanche behaviour. 
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Appendix A: Periodic Boundary Condition Treatment 

Due to the long range nature of eqn. [2j all image contributions must be taken into 
account for the correct treatment of the periodic boundary conditions along the dipolar mat 
direction. The force per unit dislocation length on a dislocation n (at position (x n , y n )) due 
to a dislocation n' (at position (x n >,y n i)) and it's images, is given by the infinite summation: 



„ ^<Vx,nUx,n' \ ^ ~ \^n' ~ KU))yyx n - yx n > - ku)) - yy n - y n >) ) 

fx ' nn ' ~ 2tt(1 - u) fe f4 {{x n - (x n , - kd)) 2 + (y n - y n ,yy ' [ } 

Due to the periodic boundary condition being one dimensional, an analytic solution to the 



above summation can be found. To do this, the summation in eqn. [ATI is re-written as 



(Ax + kd)((Ax + kdy-Ay 2 ) 
k ™™ k J^,_ x ((Ax + kd) 2 + Ay 2 ) 2 ' 1 ' 

where Ax = x n — x n > and Ay = y n — y n >, and < Ax < d/2. 

The above summation is explicitly convergent for all values of k', and may be evaluated 
analytically in terms of poly-gamma functions. This is achieved by first expressing the 



summand in eqn. A2 as an irreducible partial fraction: 



1 1 

+ -ZTl. 1 ^ + 



2 (Ax - i Ay + kd) 2 (Ax + lAy + kd) 
lAy lAy 



(A3) 



2(Ax - lAy + kd) 2 2(Ax + lAy + kd) 2 ' 

In this form, the summation to a finite k' may be straightforwardly obtained via a known 
series expansion of poly-gamma functions^, giving 

fc' 

(z + kd) d 



E 7^Y = -sk° ) (- 1 -*'^)-^ + '=' + s)l <A4) 
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and 



1 



1 

d? 



^(-l-k'+ Z -)-<f ) ^(l + k'+ Z - / 



(A5) 



(z + fed) 5 

where (f)^ m \z) is the m'th derivative of the di-Gamma function (the logarithmic derivative 
of the Gamma function^). Performing the substitutions and taking the leading order term 



as k' — > oo results in eqn. A2| evaluating to 

-d'K 



1 

2d? 



ir(Ax + tAy)\ (n(Ax-iAy) 
cot ( ; ) + cot 



d 



d 



cot 2 [ <**+^v)\ + cot2 f^ Ax - 



-in*Ay _ , rf ; , _ v d } 

For a form that is more amenable to computation the trigonometric identity, 

(a + b\ cos a + cos b 

cot = -, 

V 2 / sm a + sin b 

is used, to obtain the final result for eqn. |A2| as 

_,r sin [ rf ( cos (2^) _ cosh (2^)) +2 7rAysin(^)] 



(A6) 



(A7) 



d 2 (cos (^)- cosh 



(AS 



Eqn. A8 with the prefactor Gb x ^ n b x y /2ir(l — v) replaces eqn. ^ in eqn. when used in 



evaluating the total force of each dislocation: eqn. [3] For dislocations on the same slip 
plane, Ay = 0, and eqn. A8 reduces to the known resullP^, 



7T / TcAx 

d V d 



(A9) 



Using, eqn. |A8| with a dipolar mat periodicity length of d = 100 /mi, fig. [12] displays the 
shear-stress dependence between two dislocations of opposite Burgers vector and therefore 
on different slip planes (see fig. [T|, for different values of the dislocation density and thus 
different h. Inspection of this figure reveals the short range structure to be similar to that 
of the bare interaction given by eqn. [2] in sec. [XXJ whereas at larger distances the interaction 
correctly limits to zero at ±d/2. In this figure the chosen numerical values of the mobile 
dislocation density (where h = 1/-J~p^) will be typical of those used in the present work. 
Inspection of this figure reveals that for the highest dislocation density, the restoring stress is 
at most of the order of 1 MPa and thus small when compared to the considered values of To, 
and also typical values of the stress associated with intra-slip-plane dislocation interaction at 



sub-micron distances (see the dashed line in fig. 12, which displays the restoring shear-stress 
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FIG. 12: Plot of the restoring shear-stress as a function of separation between two dislocations of 
opposite burgers vector on different slip planes for typical dislocation densities used in the described 
simulations (coloured lines) and two dislocations of equal Burgers vector on the same slip plane 



(dashed black line). For these curves, the exact image summation form, eqn. A8 is used for a 
system with a periodicity of d = 100 fim. 



between two dislocations of the same Burgers vector on the same slip plane). To investigate 
the regime in which the two slip planes strongly interact, higher mobile dislocation densities 
would need to be considered or a different relationship between h and the dislocation density 
taken. 
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